The importance of utilizing travel history metadata for informative phylogeographical inferences: a case study of early SARS-CoV-2 introductions into Australia.

Inferring the spatiotemporal spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) via Bayesian phylogeography has been complicated by the overwhelming sampling bias present in the global genomic dataset. Previous work has demonstrated the utility of metadata in addressing this bias. Specifically, the inclusion of recent travel history of SARS-CoV-2-positive individuals into extended phylogeographical models has demonstrated increased accuracy of estimates, along with proposing alternative hypotheses that were not apparent using only genomic and geographical data. However, as the availability of comprehensive epidemiological metadata is limited, many of the current estimates rely on sequence data and basic metadata (i.e. sample date and location). As the bias within the SARS-CoV-2 sequence dataset is extensive, the degree to which we can rely on results drawn from standard phylogeographical models (i.e. discrete trait analysis) that lack integrated metadata is of great concern. This is particularly important when estimates influence and inform public health policy. We compared results generated from the same dataset, using two discrete phylogeographical models: one including travel history metadata and one without. We utilized sequences from Victoria, Australia, in this case study for two unique properties. Firstly, the high proportion of cases sequenced throughout 2020 within Victoria and the rest of Australia. Secondly, individual travel history was collected from returning travellers in Victoria during the first wave (January to May) of the coronavirus disease 2019 (COVID-19) pandemic. We found that the implementation of individual travel history was essential for the estimation of SARS-CoV-2 movement via discrete phylogeography models. Without the additional information provided by the travel history metadata, the discrete trait analysis could not be fit to the data due to numerical instability. We also suggest that during the first wave of the COVID-19 pandemic in Australia, the primary driving force behind the spread of SARS-CoV-2 was viral importation from international locations. This case study demonstrates the necessity of robust genomic datasets supplemented with epidemiological metadata for generating accurate estimates from phylogeographical models in datasets that have significant sampling bias. For future work, we recommend the collection of metadata in conjunction with genomic data. Furthermore, we highlight the risk of applying phylogeographical models to biased datasets without incorporating appropriate metadata, especially when estimates influence public health policy decision making.

Evolutionary rate of SARS-CoV-2 increases during zoonotic infection of farmed mink.

To investigate genetic signatures of adaptation to the mink host, we characterised the evolutionary rate heterogeneity in mink-associated severe acute respiratory syndrome coronaviruses (SARS-CoV-2). In 2020, the first detected anthropozoonotic spillover event of SARS-CoV-2 occurred in mink farms throughout Europe and North America. Both spill-back of mink-associated lineages into the human population and the spread into the surrounding wildlife were reported, highlighting the potential formation of a zoonotic reservoir. Our findings suggest that the evolutionary rate of SARS-CoV-2 underwent an episodic increase upon introduction into the mink host before returning to the normal range observed in humans. Furthermore, SARS-CoV-2 lineages could have circulated in the mink population for a month before detection, and during this period, evolutionary rate estimates were between 3 × 10-3 and 1.05 × 10-2 (95 per cent HPD, with a mean rate of 6.59 × 10-3) a four- to thirteen-fold increase compared to that in humans. As there is evidence for unique mutational patterns within mink-associated lineages, we explored the emergence of four mink-specific Spike protein amino acid substitutions Y453F, S1147L, F486L, and Q314K. We found that mutation Y453F emerged early in multiple mink outbreaks and that mutations F486L and Q314K may co-occur. We suggest that SARS-CoV-2 undergoes a brief, but considerable, increase in evolutionary rate in response to greater selective pressures during species jumps, which may lead to the occurrence of mink-specific mutations. These findings emphasise the necessity of ongoing surveillance of zoonotic SARS-CoV-2 infections in the future.

Genomic characterisation reveals a dominant lineage of SARS-CoV-2 in Papua New Guinea.

The coronavirus disease pandemic has highlighted the utility of pathogen genomics as a key part of comprehensive public health response to emerging infectious diseases threats, however, the ability to generate, analyse, and respond to pathogen genomic data varies around the world. Papua New Guinea (PNG), which has limited in-country capacity for genomics, has experienced significant outbreaks of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with initial genomics data indicating a large proportion of cases were from lineages that are not well defined within the current nomenclature. Through a partnership between in-country public health agencies and academic organisations, industry, and a public health genomics reference laboratory in Australia a system for routine SARS-CoV-2 genomics from PNG was established. Here we aim to characterise and describe the genomics of PNG's second wave and examine the sudden expansion of a lineage that is not well defined but very prevalent in the Western Pacific region. We generated 1797 sequences from cases in PNG and performed phylogenetic and phylodynamic analyses to examine the outbreak and characterise the circulating lineages and clusters present. Our results reveal the rapid expansion of the B.1.466.2 and related lineages within PNG, from multiple introductions into the country. We also highlight the difficulties that unstable lineage assignment causes when using genomics to assist with rapid cluster definitions.

Lymphocytic Choriomeningitis Virus Infection, Australia.

During a mouse plague in early 2021, a farmer from New South Wales, Australia, sought treatment for aseptic meningitis and was subsequently diagnosed with locally acquired lymphocytic choriomeningitis virus infection. Whole-genome sequencing identified a divergent and geographically distinct lymphocytic choriomeningitis virus strain compared with other published sequences.

The Emergence of SARS-CoV-2 Variants of Concern Is Driven by Acceleration of the Substitution Rate.

The ongoing SARS-CoV-2 pandemic has seen an unprecedented amount of rapidly generated genome data. These data have revealed the emergence of lineages with mutations associated to transmissibility and antigenicity, known as variants of concern (VOCs). A striking aspect of VOCs is that many of them involve an unusually large number of defining mutations. Current phylogenetic estimates of the substitution rate of SARS-CoV-2 suggest that its genome accrues around two mutations per month. However, VOCs can have 15 or more defining mutations and it is hypothesized that they emerged over the course of a few months, implying that they must have evolved faster for a period of time. We analyzed genome sequence data from the GISAID database to assess whether the emergence of VOCs can be attributed to changes in the substitution rate of the virus and whether this pattern can be detected at a phylogenetic level using genome data. We fit a range of molecular clock models and assessed their statistical performance. Our analyses indicate that the emergence of VOCs is driven by an episodic increase in the substitution rate of around 4-fold the background phylogenetic rate estimate that may have lasted several weeks or months. These results underscore the importance of monitoring the molecular evolution of the virus as a means of understanding the circumstances under which VOCs may emerge.

Implementing a Tailored Psychosocial Distress Screening Protocol in a Head and Neck Cancer Program.

OBJECTIVES/HYPOTHESIS: Psychosocial distress is common among patients with head and neck cancer (HNC) and is associated with poorer quality of life and clinical outcomes. Despite these risks, distress screening is not widely implemented in HNC care. In this study, we investigated the prevalence of psychosocial distress and its related factors in routine care of patients with HNC. METHODS: Data from medical records between September 2017 and March 2020 were analyzed. Psychosocial distress was measured by the National Comprehensive Cancer Network's Distress Thermometer (DT), and a modified HNC-specific problem list; depression and anxiety were assessed using the Patient Health Questionnaire-4. Descriptive statistics and logistic regression were conducted to report prevalence of distress, depression and anxiety, and factors associated with clinical distress. Implementation outcomes, including rates of referrals and follow-up for distressed patients, are also reported. RESULTS: Two hundred and eighty seven HNC patients completed the questionnaire (age 64.3 ± 14.9 years), with a mean distress score of 4.51 ± 3.35. Of those, 57% (n = 163) reported clinical distress (DT ≥ 4). Pain (odds ratio [OR] = 3.31, 95% CI = 1.75-6.26), fatigue (OR = 2.43, 95% CI = 1.1.7-5.05), anxiety (OR = 1.63, 95% CI = 1.30-2.05), and depression (OR = 1.51, 95% CI = 1.04-2.18) were significantly associated with clinical distress (P < .05). Of patients identified as distressed, 79% received same-day psychosocial evaluation. CONCLUSIONS: Clinical distress was identified in 57% of patients who completed the questionnaire, suggesting that an ultra-brief psychosocial screening protocol can be implemented in routine ambulatory oncology care, and identifies patients whose distress might otherwise go unrecognized. LEVEL OF EVIDENCE: 4 Laryngoscope, 132:1600-1608, 2022.

Metagenomic Identification of Viral Sequences in Laboratory Reagents.

Metagenomic next-generation sequencing has transformed the discovery and diagnosis of infectious disease, with the power to characterise the complete 'infectome' (bacteria, viruses, fungi, parasites) of an individual host organism. However, the identification of novel pathogens has been complicated by widespread microbial contamination in commonly used laboratory reagents. Using total RNA sequencing ("metatranscriptomics") we documented the presence of contaminant viral sequences in multiple 'blank' negative control sequencing libraries that comprise a sterile water and reagent mix. Accordingly, we identified 14 viral sequences in 7 negative control sequencing libraries. As in previous studies, several circular replication-associated protein encoding (CRESS) DNA virus-like sequences were recovered in the blank control libraries, as well as contaminating sequences from the Totiviridae, Tombusviridae and Lentiviridae families of RNA virus. These data suggest that viral contamination of common laboratory reagents is likely commonplace and can comprise a wide variety of viruses.

Genomics-informed responses in the elimination of COVID-19 in Victoria, Australia: an observational, genomic epidemiological study.

BACKGROUND: A cornerstone of Australia's ability to control COVID-19 has been effective border control with an extensive supervised quarantine programme. However, a rapid recrudescence of COVID-19 was observed in the state of Victoria in June, 2020. We aim to describe the genomic findings that located the source of this second wave and show the role of genomic epidemiology in the successful elimination of COVID-19 for a second time in Australia. METHODS: In this observational, genomic epidemiological study, we did genomic sequencing of all laboratory-confirmed cases of COVID-19 diagnosed in Victoria, Australia between Jan 25, 2020, and Jan 31, 2021. We did phylogenetic analyses, genomic cluster discovery, and integrated results with epidemiological data (detailed information on demographics, risk factors, and exposure) collected via interview by the Victorian Government Department of Health. Genomic transmission networks were used to group multiple genomic clusters when epidemiological and genomic data suggested they arose from a single importation event and diversified within Victoria. To identify transmission of emergent lineages between Victoria and other states or territories in Australia, all publicly available SARS-CoV-2 sequences uploaded before Feb 11, 2021, were obtained from the national sequence sharing programme AusTrakka, and epidemiological data were obtained from the submitting laboratories. We did phylodynamic analyses to estimate the growth rate, doubling time, and number of days from the first local infection to the collection of the first sequenced genome for the dominant local cluster, and compared our growth estimates to previously published estimates from a similar growth phase of lineage B.1.1.7 (also known as the Alpha variant) in the UK. FINDINGS: Between Jan 25, 2020, and Jan 31, 2021, there were 20 451 laboratory-confirmed cases of COVID-19 in Victoria, Australia, of which 15 431 were submitted for sequencing, and 11 711 met all quality control metrics and were included in our analysis. We identified 595 genomic clusters, with a median of five cases per cluster (IQR 2-11). Overall, samples from 11 503 (98·2%) of 11 711 cases clustered with another sample in Victoria, either within a genomic cluster or transmission network. Genomic analysis revealed that 10 426 cases, including 10 416 (98·4%) of 10 584 locally acquired cases, diagnosed during the second wave (between June and October, 2020) were derived from a single incursion from hotel quarantine, with the outbreak lineage (transmission network G, lineage D.2) rapidly detected in other Australian states and territories. Phylodynamic analyses indicated that the epidemic growth rate of the outbreak lineage in Victoria during the initial growth phase (samples collected between June 4 and July 9, 2020; 47·4 putative transmission events, per branch, per year [1/years; 95% credible interval 26·0-85·0]), was similar to that of other reported variants, such as B.1.1.7 in the UK (mean approximately 71·5 1/years). Strict interventions were implemented, and the outbreak lineage has not been detected in Australia since Oct 29, 2020. Subsequent cases represented independent international or interstate introductions, with limited local spread. INTERPRETATION: Our study highlights how rapid escalation of clonal outbreaks can occur from a single incursion. However, strict quarantine measures and decisive public health responses to emergent cases are effective, even with high epidemic growth rates. Real-time genomic surveillance can alter the way in which public health agencies view and respond to COVID-19 outbreaks. FUNDING: The Victorian Government, the National Health and Medical Research Council Australia, and the Medical Research Future Fund.

Novel hepaci- and pegi-like viruses in native Australian wildlife and non-human primates.

The Flaviviridae family of positive-sense RNA viruses contains important pathogens of humans and other animals, including Zika virus, dengue virus, and hepatitis C virus. The Flaviviridae are currently divided into four genera-Hepacivirus, Pegivirus, Pestivirus, and Flavivirus-each with a diverse host range. Members of the genus Hepacivirus are associated with an array of animal species, including humans, non-human primates, other mammalian species, as well as birds and fish, while the closely related pegiviruses have been identified in a variety of mammalian taxa, also including humans. Using a combination of total RNA and whole-genome sequencing we identified four novel hepaci-like viruses and one novel variant of a known hepacivirus in five species of Australian wildlife. The hosts infected comprised native Australian marsupials and birds, as well as a native gecko (Gehyra lauta). From these data we identified a distinct marsupial clade of hepaci-like viruses that also included an engorged Ixodes holocyclus tick collected while feeding on Australian long-nosed bandicoots (Perameles nasuta). Distinct lineages of hepaci-like viruses associated with geckos and birds were also identified. By mining the SRA database we similarly identified three new hepaci-like viruses from avian and primate hosts, as well as two novel pegi-like viruses associated with primates. The phylogenetic history of the hepaci- and pegi-like viruses as a whole, combined with co-phylogenetic analysis, provided support for virus-host co-divergence over the course of vertebrate evolution, although with frequent cross-species virus transmission. Overall, our work highlights the diversity of the Hepacivirus and Pegivirus genera as well as the uncertain phylogenetic distinction between.

The origins and genomic diversity of American Civil War Era smallpox vaccine strains.

Vaccination has transformed public health, most notably including the eradication of smallpox. Despite its profound historical importance, little is known of the origins and diversity of the viruses used in smallpox vaccination. Prior to the twentieth century, the method, source and origin of smallpox vaccinations remained unstandardised and opaque. We reconstruct and analyse viral vaccine genomes associated with smallpox vaccination from historical artefacts. Significantly, we recover viral molecules through non-destructive sampling of historical materials lacking signs of biological residues. We use the authenticated ancient genomes to reveal the evolutionary relationships of smallpox vaccination viruses within the poxviruses as a whole.

Next generation sequencing of cell free circulating tumor DNA in blood samples of recurrent and metastatic head and neck cancer patients.

BACKGROUND: Effective targeted therapy is lacking in head and neck cancer (HNC). The use of next generation sequencing (NGS) has been suggested as a way to potentially expand therapeutic options and improve outcomes. This study was performed in order to further characterize blood sample cell-free circulating tumor DNA (ctDNA) in advanced HNC patients, to determine its ability to identify actionable mutations, and to elucidate its potential role in patient management. METHODS: Retrospective analysis of 60 patients with recurrent and metastatic (R/M) HNCs who underwent molecular profiling of blood samples utilizing Guardant360, a 70-gene ctDNA NGS platform. ctDNA sequencing data was compared to tumor NGS data, when available. Best response to therapy was assessed using RECIST measures. RESULTS: The most common tumor type was oropharyngeal squamous cell carcinoma (n=21). Other cancer types included salivary gland (n=8) and thyroid (n=4). The most common mutations identified by blood analysis were TP53 (68% of patients), PIK3CA (34% of patients), NOTCH1 (20% of patients), and ARID1A (15% of patients). These findings were consistent with results from tumor sequencing data (n=30) where TP53 (48%) and PIK3CA (24%) were also the most common. Seventy-three percent (n=22) of patients had alterations identified in blood that were not present in tumor specimens. In patients with squamous cell carcinoma, 66% had an off-label option identified and 90% had a trial option identified, while 50% of patients with salivary primaries had off-label option identified and 75% had trial options identified. All patients (n=3, 100%) with thyroid primaries had off-label and clinical trial options identified. Of patients with actionable mutations, 13% (n=8) received matched targeted therapy (MTT). Three patients had stable disease (37.5%), 3 had progressive disease (37.5%), and 2 (25%) were not evaluated at the time of follow up. Of those who did not receive targeted therapy (n=21), 11 patients had stable disease (52.4%), 9 had progressive disease (42.9%), and 1 had a complete response (4.8%). CONCLUSIONS: Alterations identified by ctDNA may help inform management decisions in advanced HNC. The majority of patients had unique mutations identified on ctDNA. The role of NGS of ctDNA should be explored in future studies.

Perspectives on the Treatment of Advanced Thyroid Cancer: Approved Therapies, Resistance Mechanisms, and Future Directions.

For differentiated thyroid cancer (DTC), systemic therapy with radioactive iodine (RAI) is utilized for radiosensitive disease, while for radioiodine refractory (RAIR) disease, current standard of care is treatment with multikinase tyrosine kinase inhibitors (TKI). For BRAF-mutant DTC or anaplastic thyroid cancer (ATC), treatment with inhibitors targeting BRAF and MEK are important advances. RET-inhibitors for RET-mutated medullary thyroid cancer (MTC) recently have been FDA-approved for metastatic disease. Nevertheless, treatment of thyroid cancer resistant to current systemic therapies remains an important area of need. Resistance mechanisms are being elucidated, and novel therapies including combinations of BRAF and MEK inhibitors with RAI or other targeted therapies or TKIs combined with checkpoint inhibition are current areas of exploration.

Diversity and Evolution of Novel Invertebrate DNA Viruses Revealed by Meta-Transcriptomics.

DNA viruses comprise a wide array of genome structures and infect diverse host species. To date, most studies of DNA viruses have focused on those with the strongest disease associations. Accordingly, there has been a marked lack of sampling of DNA viruses from invertebrates. Bulk RNA sequencing has resulted in the discovery of a myriad of novel RNA viruses, and herein we used this methodology to identify actively transcribing DNA viruses in meta-transcriptomic libraries of diverse invertebrate species. Our analysis revealed high levels of phylogenetic diversity in DNA viruses, including 13 species from the Parvoviridae, Circoviridae, and Genomoviridae families of single-stranded DNA virus families, and six double-stranded DNA virus species from the Nudiviridae, Polyomaviridae, and Herpesviridae, for which few invertebrate viruses have been identified to date. By incorporating the sequence of a "blank" experimental control we also highlight the importance of reagent contamination in metagenomic studies. In sum, this work expands our knowledge of the diversity and evolution of DNA viruses and illustrates the utility of meta-transcriptomic data in identifying organisms with DNA genomes.

Comment: Characterization of Two Historic Smallpox Specimens from a Czech Museum.

The complete genome sequences of two strains of variola virus (VARV) sampled from human smallpox specimens present in the Czech National Museum, Prague, were recently determined, with one of the sequences estimated to date to the mid-19th century. Using molecular clock methods, the authors of this study go on to infer that the currently available strains of VARV share an older common ancestor, at around 1350 AD, than some recent estimates based on other archival human samples. Herein, we show that the two Czech strains exhibit anomalous branch lengths given their proposed age, and by assuming a constant rate of evolutionary change across the rest of the VARV phylogeny estimate that their true age in fact lies between 1918 and 1937. We therefore suggest that the age of the common ancestor of currently available VARV genomes most likely dates to late 16th and early 17th centuries and not ~1350 AD.

A rare case of central nervous system tuberculosis.

Intracranial abscess is an extremely rare form of central nervous system (CNS) tuberculosis (TB). We describe a case of central nervous system tuberculous abscess in absence of human immunodeficiency virus (HIV) infection. A 82-year-old Middle Eastern male from Yemen was initially brought to the emergency room due to altered mental status and acute renal failure. Cross-sectional imaging revealed multiple ring enhancing lesions located in the left cerebellum and in bilateral frontal lobe as well as in the inferior parietal lobe on the left. The patient was placed on an empiric antibiotic regimen. Preliminary testing for infectious causes was negative. Chest radiography and CT of chest showed no positive findings. He was not on any immunosuppressive medications and human immunodeficiency virus (HIV) enzyme immunoassay (EIA) test was negative. A subsequent MRI one month later showed profound worsening of the lesions with increasing vasogenic edema and newly found mass effect impinging on the fourth ventricle. Brain biopsy showed focal exudative cerebellitis and inflamed granulation tissue consistent with formation of abscesses. The diagnosis of CNS TB was finally confirmed by positive acid-fast bacilli (AFB) cultures. The patient was started on standard tuberculosis therapy but expired due to renal failure and cardiac arrest.